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ABSTRACT 

We use accurate radial velocities for 1981 member stars in 20 Galactic globular clusters, collected within our large survey aimed at 
analyzing the Na-0 anti-correlation, to study the internal kinematics of the clusters. We performed the first systematic exploration of 
the possible connections between cluster kinematics and the multiple populations phenomenon in GCs. We did not find any significant 
correlation between Na abundance and either velocity dispersion or systemic rotation. We searched for systemic rotation in the eight 
clusters of our sample that lack this analysis from previous works in the literature (NGC2808, NGC5904, NGC6171, NGC6254, 
NGC6397, NGC6388, NGC6441, and NGC6838). These clusters are found to span a wide range of rotational amplitudes from 
~0.0 km s -1 (NGC6397) to ~ 13.0 km s -1 (NGC6441). We found a significant correlation between the ratio of rotational velocity to 
central velocity dispersion (V ro </0"o) and the horizontal branch morphology parameter (B-R)/(B+R+V). The ratio V rot /o~o is found to 
correlate also with metallicity, possibly hinting at a significant role for dissipation in the process of formation of globular clusters; 
V ro , is found to correlate well with (B-R)/(B+R+V), My, (To, and [Fe/H], All these correlations strongly suggest that systemic rotation 
may be intimately linked with the processes that led to the formation of globular clusters and the stellar populations they host. 

Key words. Galaxy: globular clusters: individual: NGC104, NGC288, NGC1851, NGC1904, NGC2808, NGC3201, NGC4590, 
NGC5139, NGC5904, NGC6121, NGC6171, NGC6218, NGC6254, NGC6388, NGC6397, NGC6441, NGC6715, NGC6752, 
NGC6809, NGC6838, NGC7078, NGC7099 - Stars: abundances 



1. Introduction 

Our ideas on Galactic globular clusters (GC) have been revolu- 
tionized in the last few years by the discovery of multiple se- 
quences in their color magnitude diagrams (CMD; see Piotto 
2009; |Milone e t al. 2010 and references therein) and by evi- 
dence that the well-known star-to-star spread in the abundance 
of light elements (e.g. O, Na, Mg, Al) is not due to some mix- 
ing mechanism, but was imprinted in the gas from which cluster 
stars formed (see Rami rez & Cohen|2001[ [Gratton et al.| [2001 
2004 for discussion and references). While we lack a fully sat- 
isfying model for GC formation accounting for the variety of 
observational phenomena that are being unveiled, the notions 
that most GCs formed more than one generation of stars during 
the first few 10 8 yr of life and that the subsequent generations 
were - at least partially - polluted by chemical elements pro- 
duced by_previoiis_generatio^ accepted (see, 



e.g.|D'Antona et al |2005| |Decressin et al.|2007| 12010] | Ventura 
& D' Antona 2008 1 D'Ercole et al. 12008 [IRenzini 2008 1 Carretta 



et al. 



2010a 



Charbonnel 



[Valcarce & Catelan|2011[|Bekki|2011[|Schaerer & 



201 1| for a detailed discussion). 



Send offprint requests to: M. Bellazzini 

* Based on observations collected at ESO telescopes under programs 
072.D-507, 073.D-021 1 and 083.D-0208. 



In the previous papers of this series, the chemical signa- 
tures of this self-enrichment process were studied in detail, with 
special attention to the anti-correlation between [Na/Fe] and 
[O/Fe], using a large sample (> 2000) of stars in 19 GCs, 
for which medium-to-high resolution spectra were obtained 
with the Giraffe and UVES fiber-fed spectrographs, mounted 
at VLT@ESO (see |Carretta et al.||2006l |Gratton et al~||2007l 
|Carretta et al.||2009a|b for details, results and discussion). The 
main results of this large project have been considered in a wider 
context and correlated with other characteristics of the clusters, 
like the horizontal branch (HB) morphology, in |Carretta et al. 
( |2009c||2TjT0a] l, and |Gratton et al. | ( (20T0| >. In the present contri- 
bution we use an obvious byproduct of the spectroscopic survey, 
i.e. relatively large sets of very accurate radial velocities (RV), 
to explore the kinematics of the 19 program cluste rs (|Carretta et 
al.|2009a| Pap- VII hereafter), plus NGC1851 from |Carretta etaL 



d2011| see Table[T] below). 

These samples, which were collected for a very different sci- 
entific goal, are far from ideal to provide a full and detailed 
characterization of the cluster kinematics. The typical number 
of stars per cluster (5 150) is clearly not comparable to the most 
recent dedicated studies, based on several hundred RVs (see e.g., 
|Lane et al.|[2009| |2010a|bt |Solhma et al.||2009|) in some cases 
coupled with proper motions (van de Ven et al. 2006| van den 
|Bosch et al.|2006||McLaughlin et al.|2006| >. Furthermore, the ra- 
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dial distribution of the stars does not cover all the relevant radial 
ranges, because of limitations imposed by the fiber allocation 
process, and also because to obtain the most reliable chemical 
abundances, it is actually preferable to avoid the most crowded 
regions of the clusters where contamination from stars other than 
the intended target can become an issue. In particular most of the 
stars in our samples lie beyond the half-light radius (r/,) of the 
clusters (see Table[T] below). On the other hand, the membership 
to the parent clusters of all the stars considered in this study has 
been established not only on the basis of radial velocity cuts but 
also on very detailed chemical analysis of each individual statQ 
Moreover, our samples provide the unprecedented possibility of 
a systematic search for correlations between the chemical com- 
position and the kinematics of the cluster stars. 

The plan of the paper is the following. In Sect. [2] we briefly 
describe the global sample and assess the quality of the RV esti- 
mates. In Sect.[3]we explore the possible connection between the 
kinematics and the chemical abundance of Na, which is taken as 
the most reliable tracer of the light elements' self-enrichment at 
the basis of the Na-0 anticorrelation. In Section|4]we consider in 
more detail the RV distribution of the clusters whose kinematics 
have been previously studied only by means of integrated spec- 
tra of their center and/or with samples of RV for individual stars 
smaller than ours. In addition we search for systemic rotation in 
all the program clusters for which this analysis has not been per- 
formed in previous studies. In Sect. [5] we briefly summarize and 
discuss the main results of the analysis, including some intrigu- 
ing correlations between rotation and other cluster parameters, 
which are found here for the first time. 



2. The sample 

In the following we use the symbols V r and RV interchangeably 
to indicate heliocentric line-of-sight velocities. We make use of 
the data collected in the |Harris~| ( |1996| H96 hereafter) catalog: in 
all cases we refer to the 2010 version of this catalog, publicly 
available at www.physics.mcmaster.ca/Globular.html. 



The coordinates of the cluster centers are taken from Noyola 
|& GebhardT| ( |2006[ ), when available, and from H 96 in the other 



cases. The parameters of the King (1966 K66) or |Wilson| ( 1975 



W75) models best-fitting the surface brightness profiles of the 
clusters are taken from McLaug hlin & van der Marel | ([2005 
MvM05), specifically the core and tidal radii r c , r t , while half- 
light radii r/„ lacking in MvM05, are from H96. We adopt the 
metallicity scale by Carretta et al. (2009c). Kouwenhoven & de 
Grijs (2008) has recently demonstrated that in relatively dense 



and massive clusters like those considered here, the effects of bi- 
nary stars on the global kinematics should be minor or negligible 
(e.g. leading to an increase in the observed velocity dispersion 
S5 0.5 km s _1 , even for binary fractions as high as 100 percent, 
see their Fig. 10). Since, in addition, GCs typically have binary 
fractions ^ 20 percent ( |Sollima et al.|2007|pilone et al.|2 008), 
we consider binaries as a negligible factor for our analysi^j] 

Any detail about the measure of RV and chemical abun- 
dances, as well as target selection and membership can be found 
in the previous papers of the series ( Carretta et al.|2006 Gratton 



dense field of the Galactic bulge, distance from the cluster center was 
also taken into account when establishing membership ( Gratton et al. 



In the particularly challenging case of NGC6441, which lies in the 



2007[>. 



2 The very low number of stars showing RV differences from re- 
peated measures larg er than 2cr in the lower panel of Fig. |1| below, 
also supports this conclusion. 
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Fig. 1. Upper panel: comparison between RVs estimated from 
spectra obtained with the HR13 and HR11 grisms. Different 
symbols correspond to different clusters. The continuous line 
is drawn at a typical value for the mean difference (AV f = 
-0.5 km s~'), and the dotted lines enclose the +2<x range about 
the mean, where a typical value of cr = 0.7 km s has been 
adopted (see Table [TJ. Lower panel: comparison with the ex- 
ternal dataset by Lll for two representative clusters, NGC104 
(filled circles) and NGC6218 (empty circles). The continuous 
line marks AV r = 0.0 km s _1 , the dotted lines enclose the ±2cr 
range. The fraction of stars lying within +2<x is also reported 
[F(A<2cr)l 



|et al. |2007||Carretta et al.|2009a|b|[20TTT ). All the data used in 
this paper are publicly available in electronic form at the CDS^j 
(see the original papers). 

The program clusters are listed in Table [T] The second col- 
umn of the table reports the number of member stars having a 
valid measure of RV. The number of stars per sample ranges 
from 30 to 156 stars per cluster, with a mean of 99. Each of 
these stars also has at least a measure of iron abundance, and 
in many cases the abundance of other elements is also available 
(Na, in particular). The selected stars not only have RV that is 
fully compatible with the systemic velocity of the cluster but also 
display chemical abundance pattern identifying them as cluster 
members. The fraction of sample stars lying within the half-light 
radius of the cluster is listed in the third column of the table. 

In the large majority of clusters most stars have been 
observed using two FLAMES -GIRAFFE set-ups, the high- 
resolution gratings HR1 1 (centered on 5728 A) and HR13 (cen- 
tered on 6273 A), which were respectively chosen to measure 
the Na doublets at 5682-5688 A and 6154-6160 A, and the [O 
I] forbidden lines at 6300, 6363 A, as well as several lines of 
Fe-peak and a-elements. The spectral resolutions are R = 24200 
(for HR11) and R = 22500 (for HR13), at the center of the spec- 
tra. The fifth column of the table reports the number of stars 
having valid RV measures from spectra obtained with both set- 
ups. These stars therefore have two independent RV estimates 
and can be used to estimate the internal consistency and accu- 
racy of our radial velocities. This is especially useful since the 



|http : //cdsweb . u-strasbg . fr/ 
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pipeline used for the survey does not provide uncertainties on the 
individual RV estimates. However, it has been demonstrated by 
comparison with other independent sets of RV estimates that the 
adopted pipeline provides RV estimates with typical errors that 
are smaller than 1 k m s" 1 from spectra obtained with the same 



instrumental set-up ( Carretta et al.|201 0c p| 

The comparison between the RV estimates obtained from the 
two GIRAFFE set-ups is shown in the upper panel of Fig. [T] The 
mean difference and the standard deviation of the difference be- 
tween the two sets of estimates are reported in the 4th and 5th 
columns of Table [T] It can be appreciated that both the mean 
difference and the r.m.s. are lower than 1 km s _1 in the large 
majority of cases. The small differences in the zero-point be- 
tween the two RV sets are not a reason for concern in the present 
analysis. We transformed HR11 RVs into the HR13 system by 
applying the shifts listed in the 4th column of Table [T] then we 
obtained a single RV estimate per star by adopting RVhr\3 as 
reference value, and the corrected RVhru value when the HR13 
velocity is missing. Velocities from HR1 1 were chosen as a ref- 
erence since the associated spectra typically have higher signal- 
to-noise ratios than HR13 ones. On the other hand, we have to 
note that the HR13 - HR1 1 comparison reveals that three clusters 
have RV estimates of significantly lower accuracy than the rest 
of the sample. These are NGC4590, NGC6397, and NGC7078, 
three among the most metal-poor clusters of the sample. The 
lower quality of the RV measures for these clusters are likely 
accounted for by the extreme weakness of the spectral lines due 
to the extremely low metal content, possibly coupled with the 
relatively high temperature of the target stars. The latter effect 
depends on the actual distribution in color and magnitude of the 
targets, which can differ from cluster to cluster. For instance, 
considering the case of two very metal-poor clusters having RV 
measures of different quality, the median T e ff of the NGC6809 
sample (crA = 0.4 km s" 1 ) is 4797 K, while it is 5030 K for 
NGC7078 (<xA = 4.6 km s _1 ). Cooler atmospheres imply more 
detectable lines, keeping all the other parameters fixed, thus al- 
lowing more accurate RV estimates. 

In the lower panel of Fig. [T] we show the comparison be- 
tween the RV from the present study and those from the pub- 
lic dataset by Lane et al. ( 201 1 j Lll, hereafter) for two clus- 
ters in common between the two datasets taken as representa- 
tive of the whole range of magnitude covered by our targets, 
i.e. NGC104 and NGC6752. The zero-points agree within less 
than 0.5 km s _1 and, most important, the standard deviation 
in the RV difference is < 1.2 km s , implying typical errors 
~ 1.2/ V(2) =! 0.8 km s _1 . We note that 95%-98% of the dif- 
ferences lie within +2.0<r from the mean, which is fully consis- 
tent with the expectation for a Gaussian distribution. These re- 
sults are in excellent agreement with those presented in Carretta 
et al. (2010c I. In conclusion, internal and external comparisons 
consistently indicate that individual RV estimates in our sample 
have typical uncertainties smaller than 1 km s , so their quality 
is fully adequate for studying the internal kinematics of GCs. It 
should be kept in mind that in the cases of NGC4590, NGC6397, 
and NGC7078 the uncertainties are larger, typically 2 km s _1 , or 
even slightly larger in the case of NGC7078. 

The last column of Table Q] lists the main sources of kine- 
matic data for each cluster available in the literature, i.e. from 
the largest samples and/or the most recent and comprehensive 
studies. Eight clusters are shared with the survey by |Lane et al. 



4 See also |Pancino et al.| 1 2007} , for the accuracy of RV estimates 
obtained with FLAMES/GIRAFFE HR13 spectra by comparison with 
independent datasets. 



(2009l |2010a]bl hereafter L09, LlOa, LlOb, respectively), whose 
data have been published in Lll. These papers report compre- 
hensive kinematical analysis based on very large samples (sev- 
eral hundred member stars per cluster) of R = 10000 spectra 
(in the CaT region) obtained with the AAOmega instrument at 
the Anglo Australian Telescope, and can be considered as the 
state-of-the-art, at least from the observational point of view (see 
Lll, for discussion). A few other clusters have been studied with 



sm aller samples than in Lll but still larger than ours:|Scarpa et 
aL|( |20111 >, |C6te et al. 



(1995 , Meylan & Mayor (1991 



sembled by Pryor & Meylan (1993 PM93 hereafter). PM93 also 



provide comprehensive references for studies performed before 
1993-1994, and it remains the main source for central velocity 
dispersions (o"o) of Galactic GCs. In several cases estimates of 
<To derived from integrated spectra of the cluster center (from 
l lllingworth |1976[|Zaggia et al.|1992||Dubath, Meylan & Mayor | 
|l997[ > constitute a relevant complement to samples of RV for 
individual stars and, in some case, the only available measure 
of the cluster kinematics (NGC2808, NGC6388, NGC6441). 
NGC7078 has been the object of a very detailed and refined anal- 
ysis based on nearly two thousand RV and proper motions by |van| 
|den Bosch etaLl ( |2rj06| l. 

3. Kinematics and chemical composition 

There are only very few massive Galactic GCs whose stars dis- 
play a spread in iron (and/or calcium) abundance, either large 
(a» Cen, see |Johnson & P ilachowski 2010; Mari no et al. 201 la 



and references therein; T erzan5, |Ferraro et al.||2 009 ; Origlia et 



p]|2011| l or sma ll (M54, |Bellazzini et al.||2 008 1 Ca rretta et al 
boiObl M22, see |Da Costa et al.|2009UMarino et al |2011b|and 



references therein; NGC24 19, |Cohen et al.||2010| |Ibata et aL] 
20TT]|Cohen & Kirby|2011| NGC1851, |Carretta et al.|2011) . At 



present there is no firmly established case for a correlation be 
tween iron abundance (or metallicity, in general) and star kine- 
matics in these clusters. Interesting correlations have been re- 
ported by various authors for ai Cen (|Norris et al.|1997[ |Ferraro| 
|et al.||2002| |Soliima et aLl[2005l|Gratton et al.||20lT|7but were 



not confirmed by other studies ( |Pancino et al.|2007]|Belhni et al. 
2009), hence they remain to be further verified (see also[Bekki 



2010). Carretta et al. (201 1 ) found no difference in the kinemat- 



ics of the "metal -rich" and "metal -poor" stars in NGC1851. 

With the only exception of NGC1851, all the clusters in our 
sample are very homogeneous (internally) in iron abundance, 
with the observed star-to-star spread fully consistent with obser- 
vational errors ( |Carretta et al. 2009c Pap- VII). Therefore, they 
seem the ideal sample for a systematic search of connections be- 
tween the kinematics and the abundance of light elements known 
to display a significant spread within each cluster, and involved 
in the correlations and anti-correlations that were the main ob- 
ject of the overall program. As far as we know, this exploration 
has never been attempted before. 

All the available observational evidence on multiple popu- 
lations in GCs seems to indicate that the formation of the vari- 
ous generations of stars and the associated chemical enrichment 
must have occurred at very early epochs, approximately in the 
first few 10 8 yr of cluster life. It is therefore likely that any pos- 
sible sign of a difference in the kinematical properties between 
different generations of stars present at those times (predicted, 
for instance, in the models by Bekki 201 1 ) have been erased by 
two-body relaxation during the lifetime of GCs (see De cressirl| 
|et al.|2008) . According to the H96 catalog, the median two-body 
relaxation time of Galactic GCs is just a 1 Gyr, while the typical 
age is a 12.5 Gyr. However, the same was expected for any dif- 
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Table 1. The sample and its internal RV accuracy 



Clus. [Fe/H] Nf ol ¥{r<r h )" (My M i> 

[km s" 1 ] 



[km s" 1 ] 



N c Notes 



NGC104 

NGC288 

NGC1851 

NGC1904 

NGC2808 

NGC3201 

NGC4590 

NGC5904 

NGC6121 

NGC6171 

NGC6218 

NGC6254 

NGC6388 

NGC6397 

NGC6441 

NGC6752 

NGC6809 

NGC6838 

NGC7078 

NGC7099 



-0.76 
-1.32 
-1.16 
-1.58 
-1.18 
-1.51 
-2.27 
-1.33 
-1.18 
-1.03 
-1.33 
-1.57 
-0.45 
-1.99 
-0.44 
-1.55 
-1.93 
-0.82 
-2.33 
-2.33 



147 
110 
84 
58 
123 
149 
121 
136 
103 
33 
79 
147 
36 
144 
30 
137 
156 
39 
84 
65 



0.37 
0.39 
0.00 
0.09 
0.06 
0.48 
0.23 
0.13 
0.67 
0.00 
0.02 
0.32 
0.00 
0.35 
0.00 
0.29 
0.42 
0.49 
0.13 
0.20 



-0.38 
-0.19 

-0.81 
-0.43 
-0.09 
-1.21 
-1.78 
-0.81 
-0.07 
2.61 
-0.36 
-0.07 
-1.30 

-0.42 
-0.25 
0.19 
-0.26 
0.10 



0.24 
0.52 

1.55 
0.54 
0.87 
2.70 
0.30 
0.60 
0.30 
0.38 
0.63 
0.55 
1.59 

0.46 
0.40 
0.15 
4.57 
1.63 



84 
50 

40 
63 
56 
98 ' 
86 
63 
27 
56 
61 
23 
73" 

69 
72 
28 
58* 
30 



LlOa" 
LlOb 

INT(I76,D97), 184 stars f rom|Scarpa et al.|(|20 1 1 
INT(D97), 146 stars from |Scarpa et al.| |20 
INT(I76,Z92) 
420 stars from 



Cote et al. 



1995) 



L09 

INT(Dub97), 46 stars from |Rastorguev & Samus| ( |1991 
LlOb 

71 stars assembled from various sources by PM93^ 
LlOb 

25 stars from ] 



Rastorguev & Samus 



1991 



INT(I76) 

INT(D97), 127 stars from |Meylan & Mayor| ( |l991 
INT(I76,D97) 
LlOb 
LlOa 

> 100 stars ass embled from various sources by PM93 

{2006} 



1773 stars form 
L09 



van den Bosch et al. 



Notes. " Total number of stars with a valid RV measure. 

Fraction of sample stars lying within the cluster half-light radius r/,. 
c Number of stars having RV measure both from HR1 1 and HR13. 

d Reference to the most recent and/or largest sample of radial velocities available for the cluster, and used to study the cluster kinematics (PM93= 
|Pryor & Meylan 1993 1. INT indicates the presence of an estimate of the central velocity dispersion from integrated spectra; the associated reference 
is repo rted within brackets [176 = |Illingworth"| ( |1976| >; P89 = |Peterson et al.| < |1989| >; Z92 = |Zaggia et"al~|jl992) ; D97 = |Dubath, Meylan & Mayor| 
( |1997) l]. The INT references are reported only when large samples of individual velocities are lacking. NGC6171 andNGC6254 lack INT estimates. 
e Clusters included in the L09, LlOa and LlOb papers have been studied using samples of RV estimates several hundreds individual stars per cluster. 
We refer to that series of papers for further references. 

f Preliminary results from a sample of 107 stars have been presented by Scarpa et al. J2004f . 
* Clusters with lower quality RV estimates (see text). 



ference in the radial distribution of various generation of stars, 
while it is generally observed that Na-rich/O-poor stars are more 
centrally concentrated than Na-poor/O-rich s tarffis ee, for exam- 
ple, |D'Ercole et al.|[2008| |Decressin et"aLl[2010{ |Lardo et al 



|2011| references and discussion therein). Moreover, dedicated 
simulations are required to establish whether two-body relax- 
ation is in fact able to remove all the signs of difference in the ini- 
tial conditions over the whole radial range of the cluster ( |Lardo| 
et al.|2 01 1 1 Bekki 2010). Therefore, it is clearly worth verifying 



that any correlation between kinematic and light-element abun- 
dance is there, even if it is not expected to be observable at the 
present epoch. The main caveats for interpretating of our results 
are associated with the relatively sparse samples and by the ra- 
dial distribution of the target stars that virtually misses the clus- 
ter cores (see Table [T] also recalling that the core radii of GCs 
are always smaller than their half-light radii). 

As the best tracer of the light element spread we take the 
sodium abundance [Na/Fe]. The main reasons for this choice are 
that (a) Na abundance is much easier to measure than, for in- 
stance, O abundance, and, as a consequence, reliable Na abun- 
dance estimates are available for the large majority of the stars 
in our sample, and, (b) Na is one of the elements showing the 
widest abundance spread within each globular. 



5 Na-poor/O-rich stars are generally identified as belonging to the 
first generation of cluster stars, while Na-rich/O-poor stars are associ- 
ated to the second (or any further) generation, see e.g. Pap- VII, and 
references therein. 



In Figs. [2] and [3] we plot the velocity dispersion of cluster 
stars as a function of the [Na/Fe] abundance for all the program 
clusters. In these plots, each point corresponds to the value of the 
velocity dispersion cr estimated from the stars whose [Na/Fe] is 
enclosed within ±0.2 dex from the central values (marked by the 
plotted points); these are spaced by 0.1 dex. Only estimates ob- 
tained from > 10 stars (after clipping) are plotted. The errors 
on the velocity dispersion are computed with the jackknife boot- 
strapping technique ( Lupton|1993| l, as in Bellazzini et al. (2008 
B08 hereafter). 



In the vast majority of cases dispersion curves as a function 
of [Na/Fe] are completely flat, within the uncertainties, so the 
main conclusion that can be drawn from Figs. [2] and [3] is that, 
in general, there is no correlation between the velocity disper- 
sion and sodium abundance. There are only three cases that de- 
serve some further comment. The two metal-rich and bi-modal 
HB m orphology clusters NGC6388 and NGC6441 ( |Rich et al.| 
|1997[ ) show some marginal evidence of a drop in the veloc- 
ity dispersion with increasing Na abundance, occurring around 
[Na/Fe]~ 0.4, in line with the predictions of Bekki (2011) mod- 
els. The same seems to also occur for NGC2808 (another clus- 
ter with multimodal HB, see Dalessan dro et al.l|2011| and ref- 
erences therein), but in this case the effect is more pronounced, 
also thanks to the larger sample and the lower uncertainties in 
the cr estimates. To check the significance of the observed differ- 
ence we compared the mean RV and the velocity dispersion of 



4 



M. Bellazzini et al.: Na-0 anticorrelation and HB. IX. Kinematics of the program clusters. 
Table 2. Kinematics of Na-rich* and Na-poor* stars in selected clusters 



Clus. [Na/Fe]* r (RV Na . poor ) cr Na _ pBm N Na 

-poor 

[kms~'] [kms -1 ] [kms -1 ] [kms -1 ] 



F a , Prob." 

rip 



NGC2808 


0.4 


103.7 


10.2 


82 


99.4 


7.5 


41 


1.864 


0.008 




0.5 


103.2 


10.0 


98 


98.5 


6.6 


25 


2.269 


0.002 


NGC6388 


0.4 


84.4 


9.6 


17 


81.1 


7.2 


18 


1.732 


0.132 




0.5 


83.4 


9.3 


24 


81.1 


6.6 


11 


1.957 


0.082 


NGC6441 


0.4 


19.2 


13.7 


13 


23.6 


12.5 


16 


1.202 


0.373 




0.5 


22.2 


15.0 


17 


20.8 


10.1 


12 


2.203 


0.067 


NGC6752 


0.4 


-26.4 


5.7 


51 


-26.4 


5.0 


62 


1.302 


0.166 




0.5 


-26.4 


5.6 


60 


-26.2 


5.0 


53 


1.248 


0.202 


NGC6809 


0.4 


175.0 


4.0 


58 


174.8 


3.6 


64 


1.263 


0.184 




0.5 


174.7 


3.9 


83 


175.3 


3.5 


39 


1.254 


0.194 



Notes. * The definition of Na-rich and Na-poor: Na-rich stars have [Na/Fe]>Tr, while Na-poor stars have [Na/Fe] < [Na/Fe] 
" F p/r = H"-" di • Prob. = probability that the two samples have the same variance, according to the F-test. 
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Fig. 2. Velocity dispersion as a function of Na abundance 
for NGC104, NGC2808, NGC288, NGC1851, NGC1904, 
NGC3201, NGC4590, and NGC5904. Larger filled squares in- 
dicate cases where one star has been clipped from the sample 
because of an RV value more than 3 - <x from the mean. Only 
estimates obtained from > 10 stars (after clipping) are plotted. 
The vertical scale is different for NGC104 and NGC2808, to ac- 
commodate for the larger dispersion of these clusters. 

the Na-poor and Na-rich samples, for two different threshold^] 
separating the two subsamples, [Na/Fe] =0.4 and [Na/Fe] =0.5, 
for NGC2808, NGC6388, NGC6441, and for NGC6752 and 
NGC6809, taken as examples of clusters with flat cr vs. [Na/Fe] 
distributions. The results are reported in Table. |2) together with 

cr 2 _ 

the value of the F parameter (F p / r - X an d the proba- 

^~ Na— poor 

bility that the two subsamples are drawn from populations hav- 
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Fig. 3. The same as Fig. [ijfor the remaining clusters of our sam- 
ple. The vertical scale is different for NGC6838 and NGC6441 
to accommodate the larger dispersion of these clusters. 



ing the same velocity dispersion, according to an F test ( |Brandt 



6 This is not the same definition as adopted in Pap- VII to separate first 
and second generation stars. Here the threshold is located according to 
the observed change in the velocity dispersion. 



[T970| ). The probability is about 20 per cent for NGC6809 and 
NGC6752, independent of the adopted threshold. The possibility 
that the two subpopulations have the same cr is clearly far from 
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Fig. 4. Comparison of cumulative distributions for Na-rich (red 
continuous line) and Na-poor (blue dotted line) for NGC2808 
(upper panels) and NGC6752 (lower panels). Left panels: distri- 
butions of the radial velocity in the cluster system. Right pan- 
els: distributions of the distances from the cluster center. The 
probability that the Na-rich and Na-poor populations are drawn 
from the same parent population, according to a Kolmogorov- 
Smirnov test, is reported in the lower right corner of each panel. 
The adopted threshold between Na-poor and Na-rich stars is 
[Na/Fe]=0.5. 



excluded, for NGC6388 and NGC6441, but the probability is 
lower than ten per cent in both cases, at least for the [Na/Fe]=0.5 
threshold. Hence, the hint of a difference shown here may need 
to be explored with larger samples. On the other hand, the prob- 
ability is lower than one per cent in the case of NGC2808, inde- 
pendent of the adopted threshold. The NGC2808 sample is rich 
enough to allow for a more detailed investigation. 

3.0.1. The case of NGC2808 

An obvious possibility for the origin of the lower <x of Na- 
rich (w.r.t. Na-poor) stars in NGC2808 is that the former may 
have a radial distribution less concentrated than the latter. In 
any isotropic mass-follow-light self-gravitating equilibrium sys- 
tem, the velocity dispersion reaches its maximum value at the 
center and declines with radius. Independently of anisotropy, 
this is also observed to be the case for NGC2808 (see below). 
Although Na-rich stars are usually observed to be more cen- 
trally concentrated than Na-poor stars, still an unlucky selection 
of targets may turn out in a Na-rich sample typically lying at 
larger radii, thus implying a lower cr with respect to Na-poor. To 
check this possibility we plot in Fig.|4]the RV and radial cumu- 
lative distributions for Na-rich and Na-poor stars in NGC2808 
and NGC6752, for comparison. 

In NGC6752 the two RV distribution are indistinguishable, 
while the higher concentration of the Na-rich sample is obvious, 
and likely tracing a real difference in the radial distribution of the 
two populations (see Carretta et al.||20l"0a Lardo et al. 201 1) . 
In NGC2808 we see the opposite case: the radial distributions 



are indistinguishable, while the RV distribution of Na-rich stars 
is much steeper than for Na-poor. According to a Kolmogorv- 
Smirnov (KS) test, the probability that the two subsamples are 
from the same parent distribution of RV is just 1.1 per cent, in 
excellent agreement with the result of the F test, so the observed 
difference in cr does not come from a difference in radial distri- 
bution (the probability that the velocity distribution of Na-rich 
and Na-poor stars of NGC2808 are drawn from the same par- 
ent population decreases to 0.2% if [Na/Fe]=0.4 is adopted as 
threshold). However the upperleft hand panel of Fig.|4]provides 
an interesting hint in this sense. It can be readily appreciated that 
nearly 90 percent of the Na-rich subsample lies at negative ve- 
locity, in the reference system of the cluster (in the following 
RV or V, have always to be intended as expressed in the refer- 
ence system of the clusters, i.e. with the systemic velocity of the 
clusters subtracted). This is confirmed by looking at the mean 
velocities of the two samples, which are observed to differ by 
-4.0 km s _1 . 

In the uppermost panel of Fig. [5] it is apparent that the dis- 
tribution of RV with radius separates into two branches around 
r ~ 5 r c , as if the velocity distribution was bimodal at large radii. 
This is strongly suggestive of systemic rotation, as is indeed con- 
firmed in Sect. |4.2| Moreover, it is clear that the large majority 
of the Na-rich stars belongs to only one of the two branches. In 



Sect. 4.2 we will show that the two branches corresponds to the 



two wings of the cluster rotation curve. Likely by mere chance, 
our sample contains more stars from the half of the cluster that is 
approaching than from the half that is receding. This asymmetry 
in the distribution of targets has been exacerbated in the Na-rich 
sample, enough to produce the significant difference detected by 
the F and KS tests (that assume no selection effects, i.e. even 
sampling) but still compatible with a chance occurrence. 

In our view, there are two interesting conclusions that can 
be drawn from this case. First, it would be worth verifying this 
result with a larger and more evenly distributed sample. While 
it is likely that the observed effect is from the described rota- 
tion+sampling conspiracy, we cannot exclude that there is a real 
difference. If this were the case, it would have a striking impact 
on any scenario for the formation of this very peculiar cluster 
(P'Antona et al.|2005||Piotto et al.|2007) |Bragaglia et al.|2010[ 



Amaro-Seoane & Konstantinidis||2011|l. Second, it provides an 



example of the possible dangers connected with neglecting ro- 
tation in the analysis of velocity of stars in a GC, when the am- 
plitude of the rotation (A rot , see Sect. 4.2 1 is relatively high with 
respect to the dispersion (A rot /<To ^ 0.2). 

3.1. Rotation and Na abundance 

A distinctive prediction of the models for the formation of GCs 
with multiple populations by Bekki (2010, 201 1 ) is that even if 
a very small fraction (s=0.3 percent) of the kinetic energy of the 
gas forming the first generation (FG) of stars is in form of bulk 
rotation energy, stars in the second generation (SG) acquire sig- 
nificant rotation, much stronger than what remains in the FG. In 
particular, s values in the range 0.3 -18.0 percent correspond to 
ratios of rotational velocity to velocity dispersion (V rot /o-) from 
0.75 to 2.5 for FG stars. BekkT|( |2010l l concludes that "... the ro- 
tation of SG [stars, ndr] is due largely to dissipative accretion 
processes of AGB ejecta from FG with initially a small amount 
of rotation In Bekki's model the formation of both genera- 
tions of stars is essentially completed after a few 10 8 yr from the 
formation of the very first star, so, this possible difference in the 
rotation of different generations of cluster stars should also have 
been largely smoothed out by subsequent dynamical evolution, 
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Fig. 5. Radial velocity in the cluster system as a function of distance from the cluster center, expressed in core radii units, for less 
studied clusters. Na-rich and Na-poor stars are plotted with different symbols. The dotted line marks the tidal radius of the K66 
model that best fits the surface brightness profile, according to McLaughlin & van der Marel (2005 1. The continuous curves are the 
+3cr profiles of the same models, normalized to the central cr listed by Harris \ 1996 1, approximately representing the envelope of 
the allowed velocities for the bound members of the clusters. The dashed line marks the half-light radius, taken from Harris ( 1996 ). 



during the > 10 Gyr of a cluster life (see Decressin et al.|2 008 ). 
Nevertheless, dedicated simulations are required to check in de- 
tail whether some fossil of such difference may still be observ- 
able in GCs at the present epoch, given the initial conditions con- 
sidered in Bekki (2010). On the other hand, significant rotation 
of the cluster as a whole is a general outcome of the merging of 
two clusters ( [Making et al.|199l) , another channel that has been 
suggested as playing a role in the formation of multiple popula- 
tions GCs (Ca rretta et al.|20lT||Am aro-Seo ane & Kon stantinidis 
|201 1| and references therein). 

We searched for differences in the rotation amplitude be- 
tween Na-rich ([Na/Fe]> 0.4) and Na-poor ([Na/Fe]< 0.4) stars 
within each program cluster using the same kind of analysis 
adopted and described in Sect. 4.2 below, to look for rotation in 
the total samples. We cannot reach definitive conclusions from 
our data, since this kind of analysis would require larger sam- 
ples to allow a fully reliable study of subsamples, and we lack 
coverage of the innermost regions of the clusters, where most of 
the SG population is expected to reside in Bekki's model. In ad- 
dition, rotation amplitudes are intrinsically small in GCs (L09, 
L10a,b), and we were forced to consider only the [Na/Fe]=0.4 
threshold to avoid too poorly populated Na-rich subsamples. It 
is also worth noting that Na features are among the most promi- 
nent ones for metal-poor stars in the spectral range covered by 
our data. Therefore stars with low Na abundance (i.e. weak Na 
features) may have slightly larger uncertainties in RV than stars 
with high Na abundance. 

We simply report that we do not find any significant differ- 
ence (i.e. <: 2 km s ) in mean rotation amplitude (A rol , see be- 
low) between Na-rich and Na-poor stars for cluster having at 



least 20 stars per subsampl^] Furthermore, any marginal differ- 
ence is in the opposite sense with respect to Bekki's model ex- 
pectations; i.e., in general, Na-poor samples tend to have larger 
rotation amplitudes than Na-rich samples. 

In NGC6171 the Na-poor group displays an amplitude of 
5.5 km s _1 , but the Na-rich one only 0.5 km s . The two sub- 
samples contain only 20 and 13 stars, thus we should not trust 
this result at face value. In NGC7078 the Na-poor group displays 
an amplitude of 6.5 km s _I , and the Na-rich one just 1.0 km s -1 , 
but the two subsamples contains only 20 and 31 stars; in addi- 
tion, this is the cluster having the worst quality RV estimates in 
our survey. It is also worth noting that Na features are among 
the most prominent ones for very metal-poor stars in the spectral 
range covered by our data. Therefore stars with low Na abun- 
dance (i.e. weak Na features) may have larger uncertainties in 
RV than stars with high Na abundance, which may be relevant 
in the case of NGC7078 where the overall RV uncertainties are 
quite large. 

Finally a brief comment is needed on the clusters show- 
ing signs of Na-abundance vs. cr correlation. The samples for 
NGC6441 and NGC6388 are as small as the one for NGC6171; 
however, the Na-rich and Na-poor groups display very similar 
rotational patterns in both clusters (see below). In NGC2808, 
where the Na-rich and Na-poor groups contain 41 and 82 stars, 
respectively, Na-poor stars display a larger rotational amplitude 
than Na-rich stars {A ro , 4.5 km s~'and A rot 2.0 km s _1 , 
respectively). This is likely accounted for by the fact that the 



7 The observed differences in the position angle of the rotation axis 
maximizing the rotation amplitude (/Mo, see Sect. |4.2| l are not consid- 
ered as reliable, because of the large uncertainties associated with the 
estimate of this parameter from such small samples. 
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majority of Na-rich stars lie mainly on one of the rotation curve 



wings, as anticipated above and discussed in Sect. 4.2 



4. Dispersion and rotation for the less studied 
clusters 

There are a few clusters among those listed in Table[T]for which 
the present sample can provide a significant improvement over 
existing kinematic data and analysis. For NGC2808, NGC6388, 
and NGC6441, only estimates of o~o from integrated spectra 
are available. For NGC5904, in addition to the integrated mea- 
sure, a sample of RV for 46 individual stars was considered by 
|Rastorguev & Samus| ( |199l) , still significantly smaller than our 
sample (136 stars). Finally, NGC6254 lacks any integrated mea- 
sure, the only existing estimate of the velocity dispersion for this 
cluster coming from RVs for 25 stars from Rastorguev & Samus 
( 1991 ) in the range 0.6r e - 6.9r c , while our sample includes 147 
stars. In the following we briefly discuss the properties of the RV 
distributions in these clusters, providing new estimates of cr if 
we consider the case worth it (see below). 

In Fig. [5] we show the RV distribution in the cluster system 
as a function of projected distance from the cluster center, ex- 
pressed in units of core radius. We superposed the +3. Oct con- 
tour of the K66 model that is found by MvM05 to best fit the 
surface brightness profile of the clusters, as a reference. The po- 
sition of r/, and r, are also reportecj^] Na-rich and Na-poor stars 
are plotted with different symbols. 

At first glance, Fig. [5] reveals both (a) the (already noted) 
poor sampling of the innermost regions of the clusters provided 
by our data, and (b) the wide variety in the radial sampling from 
cluster to cluster. For example, the NGC6254 sample has a sig- 
nificant fraction of stars within r/, and barely reaches 0.5r,, while 
the innermost stars of the NGC6388 and NGC6441 lie around 
k. 2rh, and the farthest are far beyond the tidal radius. As dis- 
cussed in the introduction our samples are not best suited to a 
full characterization of the clusters kinematics. Still, in several 
cases, they are competitive with the best samples available in 
the literature, and, for the five clusters considered in this section, 
they are clearly superior to other available samples. 

The adopted K66 models appear to provide a reasonable 
description of the observed velocity distribution for NGC6254, 
NGC2808, and NGC5904. In the last two cases there are some 
stars lying beyond the +3. Oct contour of the best-fit model. This 
may be due to several trivial factors: uncertainty in the adopted 
value of cr providing the normalization constant of the model 
profile, uncertainties in the parameters of the best-fit model 
and/or limitations in the available data for the surface bright- 
ness profile (see, for example Corren ti et al. ||20TT and refer- 
ences therein), and inadequacy of the K66 model in represent- 
ing the kinematics of the clusters (e.g. rotation, anisotropy, see 
MvM05). In addition to these issues, recent detailed studies have 
pointed out that a large fraction of unbound (due to Galactic 
tides) cluster stars can contaminate the velocity dispersion pro- 
files of actual GCs well inside their Jacobi radius. For example , 
in the M — 10 4 Mq models considered by Kiipper et al. (2010 1, 
unbound stars are found to dominate the samples beyond a ra- 
dius 0.5 of the Jacobi radius. 



8 The outer edge of the K66 profiles does not exactly match the re- 
ported r, (dotted line). This is because we selected the K66 model hav- 



ing the nearest C = log(r,/r c ) to the value reported by McLaug hlin &| 
|van der Marel | {2005 ) from a grid of models sampling the range of C 
at 0.1-0.05 ste ps. For example, for NGC2808, McLaughli n & van der| 
Marel (2005 1 report C=1.56 and we adopted the model having C=1.6. 




2 4 6 

r [arcmin] 

Fig. 6. Comparison of the observed velocity dispersion profile 
for NGC5904 as derived from our data (filled circles with er- 
ror bars) and the prediction of the K66 model that best fits the 
SB profile, according to MvM05, with <x taken from H96 (thick 
continuous line). An alternative value of <x , providing signifi- 
cantly better representation of the observed profile, is proposed 
(<xo = 7.5 km s , from our own fit, dashed line). 



The NGC6388 and NGC6441 samples are too sparse to let 
us draw firm conclusions. At face value they show a pretty flat 
dispersion curve over the considered radial range. This, along 
with the significant rotation (see below) and the presence of sev- 
eral members beyond the K66 tidal radius, may be suggestive of 
strong effects from Galactic tides ( Mu noz et al.||2008| [Sollima 
et al. 2009; Kupper ^t al.|2010"l and references therein; see also 
the discussion in Sect. [5] below). The large number of members 
beyond the tidal radius seen in NGC6441 may simply be due 
to the inadequacy of K66 to describe the outer regions of clus- 
ters that is discussed in detail by MvM05. These authors found 
that in most cases spherical and isotropic W75 models (more 
extended than K66 ones) provide the best overall fit of cluster 
observed light profiles. This is actually the case for NGC6441, 
while for NGC6388, K66 is better. The tidal radii for the best-fit 
W75 models are 4.2 times and 6.1 times larger than their K66 
counterparts, for NGC6388 and NGC6441, respectively. All the 
stars in our samples would thus lie well within the W75 tidal 
radius of the clusters. 



4.1. Newer,) estimate for NGC5904 

In the comparisons with other self-gravitating systems, the kine- 
matic properties of a GC are usually summarized by a single 



number, i.e. the central velocity dispersion cr ( |DjorgovsE 1995 
Bellazzini et al.||T998] >. This is the number usually reported in 



GC catalogs (PM93, H96), so it may be useful to provide re- 
finements of such estimates or to confirm existing estimates on 
significantly sounder basis, when possible. 

We used our data to produce velocity dispersion curves for 
the less studied clusters listed above, in the same way as de- 
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Fig. 7. Rotation in NGC2808, NGC6388, and NGC6441. The 
plots display the difference between the mean velocities on each 
side of a cluster with respect to a line passing through the clus- 
ter center with a position angle PA (measured from north to 
east, north=0°, east=90°), as a function of the adopted PA. The 
adopted line coincides with the projected rotation axis. The con- 
tinuous line is the sine law that best fits the observed pattern. 
The vertical scale used for NGC6441 is larger, to accommo- 
date the strong rotation signal. It is also worth noting that even 
limiting the NGC6441 sample to the 21 stars lying within the 
tidal radius (see Fig. [5} the rotation pattern is very similar, with 
A rnt =i 1 1 km s'and PA n -225°. 



scribed in B08, using jackknife resampling to compute uncer- 
tainties. We compared the observed curves with the predictions 
of the K66 models plotted in Fig. [5] normalized with the ctq 
listed in H96, to check that our new results are compatible with 
these values. 

The agreement is acceptable or satisfying for NGC6388, 
NGC2808, and NGC6254, so we confirm the cr estimate re- 
ported in the literature for these clusters. From the sparse sample 
of NGC6441 we were able to obtain only one point of the disper- 
sion curve within the tidal radius of the best-fitting K66 model; 
this point is compatible with the listed value of cr . Outside of r t 
the velocity dispersion remains fairly large (cr ~ 10 km s~'). 

A direct comparison of the observed velocity dispersion pro- 
files with the predictions of the K66 models is plotted in Fig. [6] 
for NGC5904, the only case where we found a mismatch be- 
tween the current normalization of the velocity profile of the 
model and our data. The value listed by H96 {ctq = 5.5 + 
0.4 km s , see Tableland PM93 for the original data sources) 
is clearly too low to match the observed dispersion curve. The 
adopted <x () = 7.5 ±1.0 km s instead provides a satisfying fit, 
and we propose this value as a new, more robust estimate of the 
central velocity dispersion for this cluster. It has to be recalled 
that the accuracy of this estimate also depends on how good and 
appropriate the adopted K66 models are (see MvM05). 



4.2. Rotation 

Globular clusters are generally considered as classical examples 
of pressure-supported systems, in which rotation is nonexisting 
or negligible. The classical case of rotating GC is u> Cen, which 
has been considered as very peculiar in many aspects for a long 
time, and it displays a large-amplitude (~ 8 km s _1 ) rotation 



pattern coupled with high isophotal ellipticity (Merritt, Meylan 



& Mayor 1997 van de Ven et al. 2006 Sollima et al. 2009). 



However, recent studies based on large samples are revealing 
that rotation of typical amplitude 2-4 km s" 1 is far from un- 
common in Galactic GCs (see Cote et aT||1995| |Anderson &| 
|King||2003||van den Bosch et al.||2006| L09, L10a,b). Rotation 
amplitudes detected from RV samples are just lower limits to 
the true 3D amplitude, because of projection onto the plane of 
the sky. This may support the suggestion by MvM05 that W75 
models may be more appropriate for describing actual GCs, as 
these models can account for rotation (see Solli ma et al.[ [2009 
for a recent application; however, Wilson's models adopted by 
MvM05 do not include rotation). 

We used our sample to search for rotation in all the clusters 
for which this characteristic has not been considered in previous 
studie^] In addition to these, we also analyzed NGC1851 and 
NGC1904, whose rotation has been considered by jScarpa et al.| 
(2011), but only separated in two radial bins, while we aim at 
global properties to make proper comparisons among clusters. 



Our results are in good agreement with those by |Scarpa et al. 
( |201 To get a rotation amplitude value that is homogeneous 
with the other clusters ,we also studied the case of NGC7078, 
whose kinematics have been analyzed in much greater detail by 
|van den Bosch et al.| ( [200 6). We limit our analysis to rotation 
patterns that can be detected from radial velocities, since the data 
on rotation in the plane of the sky are only available for a handful 
of clusters (|And erson & King 2003] |van de Ven et al.|2006| [van] 
|den B osch et al. 2006|l. 

To look for rotation we used the same method as was adopted 
by |C6te et aL[ ( [1995) , [Pancino et al.| p007] ), L09, and L10a,b. 
For any given cluster, we divided the sample in two, choosing 
stars lying on the opposite side of a line passing from the cluster 
center, and we computed the difference in mean RV between 
the two subsamples. Then we rotate the dividing line by a fixed 
amount (10° or 20°, depending on the number of stars of the 
considered sample) and repeated the computation. Recording the 
difference in the mean RV for each position angle (PA) of the 
dividing line we can plot the first as a function of the second: 
a coherent sinusoidal pattern is a clear sign of rotation. In the 
adopted approach PA is defined to increase anti-clockwise in the 
plane of the sky from north (PA=0°) toward east (PA=90°). We 
fit the observed patterns with the simple sine law: 



A{V r ) = A rot sin(PA + O) 



(1) 



where <E> = 270° - PAq, and PAq is the position angle of the 
dividing line corresponding to the maximum rotation ampli- 
tude (in degrees), coinciding with the rotation axis, and A rot is 
two times the actual mean amplitude (in km s _1 ; see LlOa). 
However, it must be noted that, as the (RV) difference is aver- 
aged over the full range of radius covered by the sample, A rot /2 
should be an underestimate of the maximum rotation amplitude, 
which should be the most informative global rotation parameter 
(Meylan & Heggie 1997), since the amplitude generally varies 



9 For the remaining clusters we note that the results of our rotation 
analysis are in good agreement with those reported in the literature and 
based on larger samples, in most cases from LlOb. 
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Fig.8. The same as Fig. Fj but for NGC1851, NGC1904, 
NGC5904, NGC6254, NGC6171, NGC6397, NGC6838, and 
NGC7078. The kinematics of some of these clusters has been 
studied with samples larger than ours, still they lack a study of 
rotation (NGC6171, NGC6397, NGC6838) or an estimate of the 
rotation amplitude homogeneous with those obtained here. 



significantly with distance from the cluster center (see, for exam- 
ple, the case of a> Cen Solli maet al.|2009") . By inspecting some 
well-populated rotation curves of GCs, we noted that A rot is a 



reasonabl proxy for the actual maximum amplitude (see |Pancrno 
et al.|2007] >; in the following we use this value as an estimate of 
the rotation amplitudes and, consequently multiply by two the 
values reported by LI 0b, to make the two sets of measures fully 
homogeneous. We explored the reliability of the derived param- 
eters by repeating the analysis on subsamples and by compar- 



Fig.9. Rotation curves for NGC2808, NGC6388, and 
NGC6441. Filled and open symbols are adopted for Na- 
rich and Na-poor stars, respectively. Left panels: RV in the 
system of the cluster as function of distance from the center 
projected onto the axis perpendicular to the best-fit rotation axis 
found in Figs. 7] and [8] Stars inside and outside the tidal radius 
are plotted with different symbols for NGC6441 and NGC6388. 
All the stars in the NGC2808 sample lie within the tidal radius 
on the cluster. Right panels: comparison of the cumulative 
RV distributions of stars having X(PAq) > 0.0 (continuous 
lines) and X(PA ) < 0.0 (dotted lines). The probability that the 
two distributions are drawn from the same parent population 
(according to a KS test) is reported in each panel. We show 
rotation curves only for the two clusters having Pks < 10% plus 
NGC2808. 



ing it with external samples (e.g. for clusters in common with 
Lll), and we concluded that the estimates of A ro , are fairly ro- 
bust (we assume a typical lcr uncertainty of 0.5 km s for the 
largest samples, similar to LlOb), while PA is quite sensitive to 
the adopted sample and should be considered as uncertain at the 
±30° level in the best cases. 

The results of the analysis are displayed in Figs. [7] and [8] 
It is evident that the considered clusters cover a wide range of 
rotation properties, from no rotation (NGC6397) to an ampli- 
tude of more than 10 km s (NGC6441). It is interesting to note 
that three among the most peculiar clusters in terms of multiple 
populations, NGC2808, NGC6388, and NGC6441, display par- 
ticularly strong rotation patterns, possibly hinting at a similarity 
with some galactic nuclei ( Sefh et aT] |2008). Other clusters for 
which we detect for the first time a significant amplitude of ro- 
tation are NGC5904 and NGC6171. The sig nificant amplitude 
found for NGC7078 agrees with the results by |van den Bosch et] 
|aLl( |2006l >. 

In Fig. [9] we show the rotation curve for the cases of 
NGC2808, NGC6388, and NGC6441, i.e. the distribution of RV 
as a function of the projection of the star position on the axis 
perpendicular to the axis of rotation [X(PAq)]. In the righthand 
panels, the RV distribution of stars lying on opposite sides with 
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respect to the rotation axis are compared. When there is no ro- 
tation, the two distributions should be indistinguishable, while 
with significant rotation a shift should be apparent. The degree 
to which the two distributions differ also depends on the ratio 
between rotation and velocity dispersion and on the actual shape 
of the rotation curve. In the considered cases the comparison is 
not fully conclusive, given the small samples under considera- 
tion. A KS test states that it is relatively unlikely (but clearly not 
impossible) that the observed patterns may emerge by chance 
from nonrotating systems. For NGC2808 we had the opportu- 
nity to fully confirm the result presented here with a much larger, 
still unpublished sample of RV estimates (more than 800 stars, 
Bragaglia & Carretta, private communication). The rotation pat- 
tern observed in NGC1851 has also been independently con- 
firmed with an unpublished sample of ~ 100 HB stars observed 
within another program (Gratton & Lucatello, private communi- 
cation). 

Figure [9] deserves a couple of further comments: 

- In the upper lefthand panel of Fig. [9] the subtle selection ef- 
fect of the uneven spatial sampling in NGC2808 can be ap- 
preciated. The effect of rotation (that favors a certain range 
of velocities for stars in a given wing of the rotation curve) 
coupled with low-number statistics may well be at the origin 
of the correlation between Na-abundance and velocity dis- 
persion discussed above. In fact, in each wing the velocity 
dispersion is lower than in the two wings considered as a 
whole. 

- It is interesting to note that in NGC6441 stars within and out- 
side the tidal radius seem to trace the same rotation pattern. 
The kind of solid body rotation extending to large distances 
from the cluster center that is hinted at in Fig.|9]may be sug- 
gestive of a velocity gradient induced by tidal stripping (see 
IMufioz et al.|20081|Bellazzini et al.|2008]|Sollima et al.|2009| 
discussion and references therein). 

We have to stress again that the results presented in this sec- 
tion may also suffer from biases associated with the dimension 
of the samples and with the radial distribution of sample stars. 
However, most of the stars at relatively large distances from the 
cluster centers may not be a strong concern for detecting rota- 
tion, since it is expected that the highest rotation amplitudes are 
reached relatively far from the center (Merritt, Meylan & Mayor 
1997 Scarpa et al. 201 1[ ). Moreover, the adopted technique is 
not expected to be particularly sensitive, since any rotation curve 
is strongly smoothed out in the simple average velocity of two 
halves of the clusters: as a result, the detected signals should be 
considered as quite robust and homogeneous. 

5. Summary and discussion 

We used the radial velocity estimates obtained from a large spec- 
troscopic survey (aimed at the study of the Na-O anticorrelation 
in GCs) to study the kinematics of the surveyed clusters, with 
the main aim of looking for correlations between kinematics and 
chemical abundance of elements involved in the anticorrelation 
issue. With all the limitations associated with the relatively small 
samples and the uneven radial distribution of target stars within 
each cluster, the data considered here allow the exploration of 
this interesting issue for the first time. By internal and external 
comparisons we have demonstrated that the typical uncertainty 
on individual RV estimates is smaller than 1 km s _1 . 

In general, we did not find any evidence of correlation of Na 
abundance with velocity dispersion or rotation within each clus- 
ter. A possible decrease in velocity dispersion for [Na/Fe]^ 0.4 



for the peculiar clusters NGC6388 and NGC6441, in broad 



agreement with the expectations from the models by Bekki 
( |2011| l, clearly requires larger datasets to be confirmed or dis- 
confirmed. A statistically significant decrease in velocity disper- 
sion for [Na/Fe]<; 0.4 forNGC2808 is found to likely come from 
an uneven spatial distribution of the target stars in this cluster. 

We investigated in slightly deeper detail the cases of the clus- 
ters whose kinematics were previously studied in the literature 
with smaller samples than ours. We confirm (or provide support 
for) the central velocity dispersion estimates reported in the lit- 
erature for NGC6388, NGC2808, and NGC6254. We provide a 
revised estimate of the central velocity dispersion for NGC5904. 

Finally we looked for systemic rotation in the clusters that 
were lacking such analysis in the literature (or had been studied 
with different methods than that adopted here). We found that the 
considered clusters cover a remarkable range of rotational ampli- 
tudes, from ~ km s~'to ~ 13 km s _1 . NGC5904, NGC6171, 
NGC2808, NGC6388, and NGC6441 show especially strong ro- 
tation patterns (A rot ^ 2.5 km s _1 ). 

In Table [3] we list homogeneous o~q and A ro , estimates col- 
lected from this work and other sources in the literature (mainly 
LI 0b), together with other relevant cluster parameters from var- 
ious sources (mainly from H96). We added to our original sam- 
ple co Cen (NGC5139), M54 (NGC6715), M22 (NGC6656), and 
M53 (NGC5024) for which we had all (or most of) the relevant 
data in the same scales adopted here from other sources (LI 0b; 
B08; Carretta et al.|20i0b"l|Pancino et al.|2007") . It is interesting 
to check that the rotation properties of the clusters do correlate 
with some other characteristics, as this is the largest collection of 
rotational amplitudes for Galactic GCs currently available (see 
Mey ian & Heggie|[T997l for a summary of the data available 
circa one decade ago). Before going on with the analysis it may 
be worth recalling that projection effects should be at work, and 
the listed A ro , are only estimates of the amplitude of the pro- 
jected rotation pattern, thus lower limits of the actual rotation 
amplitude. However, in Appendix\A\we show that A ro , is a rea- 
sonable proxy for the true amplitude, in a statistical sense (see 
also Chandrasekhar & Munch 1949). 

As the most informative parameter, we use the ratio between 
the rotation velocity V rol , here represented by A rnt , and the cen- 
tral velocity dispersion ctq. In Fig. 
A r< ,t/o-o vs. the inter quartile range ( 



10 we show scatter plots of 
a parameter defined by 



QR) 

|Carretta| ( |2006| l to measure the extension of the Na-O distribu- 
tion in clusters, the cluster ellipticity, the absolute integrated V 
magnitude {My), the iron abundance (metallici tv, [Fe/HD , the 



r [To 

~+v\ 

f (foe 



(Lee 



1990), the 



classical HB morphology parameter 

logarithm of the central luminosity density (togpo), the absolute 
value of the distance from the Galactic plane (|Z|), and the dis- 
tance from the Galactic center (Rcc)- 

The ratio A rof /cr does not display any clear correlation with 
My, ellipticity, logpo, |Z|, and Rcc- The lack of a connection with 
ellipticity may appear somehow surprising. However, it must be 
recalled that the two quantities may refer to different regions 
of the clusters, since it is hard to reliably estimate ellipticity at 
large distances from the cluster core. On the other hand, a hint 
of a trend with IQR is perceived. In general clusters with higher 
Arotlo'd ratios have smaller IQR. Four of the five clusters lying 
on a nearly parallel sequence above the bulk of the other objects 
(at larger IQR for a given A ro ,/cr value) are among the most pe- 
culiar clusters, two of them also displaying a spread in the iron 
abundance. According to a two-tailed Student's test, the proba- 



10 Where B,R, and V are the number of HB stars lying to the blue, to 
the red, and within the RR Lyrae instability strip, respectively. 
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Table 3. Cluster parameters 



Clus. 


0"o 
km s 


km s 


km s 


5* 

km s 


[Fe/H] 


IQR|Na/OJ 


B-R 
B+R+V 


M v 


ell 


logpo 

i 

Lqpc 


Rgc 
kpc 


kpc 


Notes 


NGC104 


9.6 


0.6 


4.4 


0.4 


-0.76 


0.472 


-0.99 


-9.42 


0.09 


4.88 


7.4 


-3.1 


LlOb 


NGC288 


2.7 


0.8 


0.5 


0.3 


-1.32 


0.776 


0.98 


-6.75 


— 


1.78 


12.0 


-8.9 


LlOb 


NGC1851 


10.4 


0.5 


1.6 


0.5 


-1.16 


0.693 


-0.32 


-8.33 


0.05 


5.09 


16.6 


-6.9 


H96+t.w. a 


NGC1904 


5.3 


0.4 


0.6 


0.5 


-1.58 


0.759 


0.89 


-7.86 


0.01 


4.08 


18.8 


-6.3 


H96+t.w. a 


NGC2808 


13.4 


1.2 


3.3 


0.5 


-1.18 


0.999 


-0.49 


-9.39 


0.12 


4.66 


11.1 


-1.9 


H96+t.w. 


NGC3201 


4.5 


0.5 


1.2 


0.3 


-1.51 


0.634 


0.08 


-7.45 


0.12 


2.71 


8.8 


0.7 


C95 


NGC4590 


2.4 


0.9 


1.2 


0.4 


-2.27 


0.372 


0.17 


-7.37 


0.05 


2.57 


10.2 


6.0 


LlOb 


NGC5024 


4.4 


0.9 


0.0 


0.5 


-2.06 


— 


0.81 


-8.71 


0.01 


3.07 


18.4 


17.6 


LlOb 


NGC5139 


19.0 


1.0 


6.0 


1.0 


-1.64 


0.930 


— 


-10.26 


0.17 


3.15 


6.4 


1.3 


VV06+P07 


NGC5904 


7.5 


1.0 


2.6 


0.5 


-1.33 


0.741 


0.31 


-8.81 


0.14 


3.88 


6.2 


5.5 


t.w. 


NGC6121 


3.9 


0.7 


1.8 


0.2 


-1.18 


0.373 


-0.06 


-7.19 


0.00 


3.64 


5.9 


0.6 


LlOb 


NGC6171 


4.1 


0.3 


2.9 


1.0 


-1.03 


0.522 


-0.73 


-7.12 


0.02 


3.08 


3.3 


2.5 


H96+t.w. c 


NGC6218 


4.7 


0.9 


0.3 


0.2 


-1.33 


0.863 


0.97 


-7.31 


0.04 


3.23 


4.5 


2.1 


LlOb 


NGC6254 


6.6 


0.8 


0.4 


0.5 


-1.57 


0.565 


0.98 


-7.48 


0.00 


3.54 


4.6 


1.7 


H96+t.w. 


NGC6388 


18.9 


0.8 


3.9 


1.0 


-0.45 


0.795 


-0.65 


-9.41 


0.01 


5.37 


3.1 


-1.2 


H96+t.w. 


NGC6397 


4.5 


0.6 


0.2 


0.5 


-1.99 


0.274 


0.98 


-6.64 


0.07 


5.76 


6.0 


-0.5 


MM91+t.w. 


NGC6441 


18.0 


0.2 


12.9 


2.0 


-0.44 


0.660 


-0.76 


-9.63 


0.02 


5.26 


3.9 


-1.0 


H96+t.w. 


NGC6656 


6.8 


0.6 


1.5 


0.4 


-1.70 




0.91 


-8.50 


0.14 


3.63 


4.9 


-0.4 


LlOb 


NGC6715 


16.4 C 


0.4 


2.0 


0.5 


-1.56 


1.169 


0.54 


-9.98 


0.06 


4.69 


18.9 


-6.5 


109 + B08 


NGC6752 


5.7 


0.7 


0.0 


0.0 


-1.55 


0.772 


1.00 


-7.73 


0.04 


5.04 


5.2 


-1.7 


LlOb 


NGC6809 


2.7 


0.5 


0.5 


0.2 


-1.93 


0.725 


0.87 


-7.57 


0.02 


2.22 


3.9 


-2.1 


LlOb 


NGC6838 


2.3 


0.2 


1.3 


0.5 


-0.82 


0.257 


-1.00 


-5.61 


0.00 


2.83 


6.7 


-0.3 


H96+t.w. 


NGC7078 


13.5 


0.9 


3.8 


0.5 


-2.33 


0.501 


0.67 


-9.19 


0.05 


5.05 


10.4 


-4.8 


H96+t.w. 


NGC7099 


5.0 


0.9 


0.0 


0.0 


-2.33 


0.607 


0.89 


-7.45 


0.01 


5.01 


7.1 


-5.9 


LlOb 



Notes. The references in the last column indicate the source of the velocity dispersions and of the rotation amplitudes (t.w.= this work; H96= 
2010 version of me [HamT|([T996ll cat alog; B08 =|Bellazzini et al |<[200"8} ; I09= |Ibata et aT7| <[2009 ); P07= |Pancino et al.|p057) ; MM91= |Meylan &| 
|Mayor| ( [T99T) ; C95 HC6te et al.|H1995) ; VV06 = |van de Ven et al.H2006} ; for the other acronyms see Tab.|l|l. T he A rot values from LlOb have been 
multip ied by 2 to report them on the same scale as ours. The metallicity values are from Carretta et al. ( 2009c|, except for NGC 1851 (from |Carretta| 
et al.|2011> and NGC 6715=M54 (from |Carretta et al.|2010c|. IQR[Na/Q l values are taken from lCarretta et al.|<[2010a) except for NGC 6715 (from 
Carretta et al.l2010b|l and NGC5139=a> Cen (from Carretta et al. 201 lb. M v and ell values are from H96. 7^1^77 values are taken from Mackey & 



van den Bergh (2005 1 except for NGC6388 and NGC6441; tor these cluster we computed updated values (properly accounting for the extended 
blue HB tail) from the data reported in Table 3 of Busso et al. ( 2 007 ) 



a In agreement with the results by |Scarpa et al. 1 201 1). " In agreem ent with the results shown in Scarpa et al. 1 2004| >. c Excluding the inner cusp 
possibly associated with a central Black Hole (see |Ibata et al. 12009] for details) 



bility that a Spearman rank correlation coefficient equal or larger 
(in absolute value) than the observed one (s = -0.353) is pro- 
duced by chance from uncorrected quantities is P, =12.6 per 
cent, so the correlation cannot be considered as statistically sig- 
nificant ( |Press et"aL]|1992| see Table |4] for s and P, for all the 
considered parameters). For this reason we refrain from any in- 
terpretation of this plot, which can be considered as a possible 
starting point for further analysis. 

An interesting trend of A rot /<To is with the cluster metallicity: 
clusters with higher metallicity have greater fractions of ordered 
motion with respect to pressure support. In spite of the three 
outliers, the correlation is found to be remarkably significant by 
the same kind of Student test as used above, P, =1.5 per cent. 
Since larger metal content in a gas implies higher efficiency in 
energy dissipation by atomic transitions, the observed correla- 
tion, if confirmed, may hint at a significant role of dissipation in 
the process of cluster formation (see Bekki 2010|[20TT l. 

The most significant correlation of A, o ,/<x , however, is with 



the HB morphology ( 



B-R 



): the relevance of ordered motions 



with respect to pressure is stronger for clusters with redder HBs. 
The probability that the observed correlation arose by chance is 
negligible (P, = 2xl(T 5 ). It is well known that an( j [Fe/H] 

are correlated (since metallicity is the first parameter determin- 
ing the HB morphology, see Lee 1990| [Fusi Pecci et al.| [l993; 



Grattonetal. 2010 and references therein), so it is likely that the 
correlation of one of the two parameters with A rof /cr is a sec- 
ondary effect of the correlation with the other. We cannot provide 
a firm identification of the primary correlation; however, we note 
that for our sample the correlation of A ro/ /<xo with B +Z+ v (s=- 
0.769) is much stronger than both that of A ro ,/<To with [Fe/H] 
(s=0.492) and of w i tn [Fe/H] (s=0.661). As a result, the 

connection seems more direct with the HB morphology than 
with the metallicitypl It is interesting to note that B ^l v shows 
some degree of correlation {\s\ ~ 0.3; not shown here) with other 
HB morphology parameters considered in |Carretta et al. (2010a 



like e.g., log r™^ Hg ) and/or estimates of Helium abundance (e.g 



Gratton et al. |2010| , in the sense that clusters with larger 
logT'™f HB and Y max have lower A ro ,/o"o. These parameters, in 
turn, positively correlate with IQR that, as seen above, displays 
an anticorrelation with A rot /cro- In summary, the overall observa- 



11 A possible connection between cluster systemic rotation and 
HB morphology was already suggested more than 25 years ago by 



|Buonanno et al.| J1985) , inspired by a previous suggestion by |Norris| 
4 1 983[ >. However, the chain of processes envisaged by these authors was 
expected to produce bluer HB morphologies for greater degrees of sys- 
temic rotation, i.e. a correlation in the opposite sense with respect to 
what is found here. 
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Fig. 10. The ratio between the amplitude of the rotation A ro , to the central velocity dispersion cr is plotted versus various other 
parameters. The clusters that appear to lie at odd positions with respect to the main trends are labeled. In the rightmost panel of the 
lower row, clusters having [Fe/H]> —1.0 are plotted as filled circles. 



tional scenario strongly suggests a significant link between the 
systemic rotation of the clusters and the processes that lead to 
the formation and evolution of the different generations of clus- 
ter stars, whose effects are likely imprinted (and most evident) 



in the current HB morphology (|Fusi Pecci et al. 1993 D'Antona 



|et al.|2005)|Gratton et al. |2010| and references therein) 

It should also be noted that the orbits of the most metal- 
rich (and HB-red) clusters are likely confined within the Galactic 
bulge, so they may be subject to stronger tides, possibly produc- 
ing velocity gradients that can mimic rotation (see e.g., Munoz et 



|al.|2008[ ). However, in this case one would expect to see a better 
correlation between cluster position and A rot than with A rot /cro, 
since the amplitude of the gradient should mainly depend on the 
orbital parameters. The Spearman rank correlation coefficients 
listed in Table [5] show that this is not the case; the correlation is 
poor also between A ror and the orbital parameters of the clusters 
(Rper and Z max from |Dinescu et al.|1999 Casetti-Dinescu et al. 
[20071 see Table [5). 

On the other hand, A rot displays significant correlations with 
v , My, cr , and [Fe/H] (P, < 1.0 percent, in all cases; see 
[TT] ). Also in this case, the strongest and most significant 
correlation is, by far, with HB morphology, therefore A rot seems 
strictly tied with intrinsic cluster parameters, not with their or- 
bit or environment. Note that My and cr are also correlated be- 
tween each other, both approximately tracing the cluster ranking 
in mass, but do not appear to correlate with [Fe/H] and/or B f fi f v ■ 
The only outlier of the correlations of A rot with My and ctq is 
NGC6441, which has a much larger A rot for its My or <x than 
do the other clusters of the sample. This may support the idea 
of a tidal origin for the strong velocity gradient observed in this 
specific cluster, as already suggested in Sect. 4.2 For the other 



B-R 
B+R+. 

Fig, 



clusters rotation is more likely of internal origin, possibly linked 
with the phase of cluster formation, with tides as a possible 
source of noise in the correlations. Removing NGC6441 from 
the sample has a negligible effect on the Spearman coefficients 
and on the statistical significance of the correlations, as mea- 



sured by P,. We found two linear combinations of parameters 
displaying correlations with A rot that are significantly stronger 
than all those with single parameters, LCI - gfsfy + 0.47 My 
and LC2 = My - 1.73 [Fe/H]. The significance of both correla- 
tions is very high (P, < 10~ 5 ). It is interesting to note that My 
(likely as a proxy of the total cluster mass) has been found to play 
a role in bi-variate correlations, both with some HB morphology 
parameters ( |Fusi Pecci et al. 1993 Recio-Blanco et al. 2006) 
and with parameters related to the anticorrelation phenomenon 
(Pap-VII, |Carretta et al.|2009b> . 

As typical of most exploratory studies, our analysis does not 
provide firm conclusions, but calls for further investigations and 
perhaps opens some new windows on the research on globular 
clusters. It is clear that a follow up of the possible chemistry- 
kinematics connection and of the rotation patterns detected here 
is necessary for NGC2808 and, above all, for NGC6388 and 
NGC6441. The possible correlation between A, o ,/cr and IQR 
also deserves further study. This requires a considerable obser- 
vational effort, e.g, to double the number of GCs with relatively 
large sample of Na abundances and RV estimates. The corre- 
lations between A rot /o"o, J£j? Y and [Fe/H] should be easier to 
verify and/or extend, since only large RV samples are required, 
with no need of detailed chemical abundance analysis of indi- 
vidual stars. The same is true for the correlations between A rot 
and 2jfify> My, cr , and [Fe/H]. These, together with those of 
A r ot/cro with HB morphology and metallicity, are especially in- 
teresting since they are quite significant, from a statistical point 
of view, and far from trivial, since they indicate that rotation may 
be intimately linked with the internal physics of the clusters, and, 
possibly, with their origin. 
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Table 5. Spearman rank correlation coefficients for A r , 



Parameter 


Spearman s" 




N c , 

cuts 


Pf 


(B-R)/(B+R+V) 


-0.775 


5.614 


23 


lxl0~ 5 


My 


-0.611 


3.623 


24 


0.001 


(To 


0.604 


3.558 


24 


0.002 


[Fe/H] 


0.531 


2.941 


24 


0.007 


ell 


0.358 


1.755 


23 


0.092 


logpo 


0.245 


1.186 


24 


0.248 


per 


-0.244 


1.097 


21 


0.286 


IQR 


0.117 


0.526 


22 


0.605 


Rac 


-0.109 


0.512 


24 


0.614 


7* 
max 


-0.106 


0.467 


21 


0.646 


z 


-0.048 


0.224 


24 


0.825 


LCI" 


-0.888 


8.865 


23 


2x1 0~ 8 


LC2 b 


-0.811 


6.509 


24 


IxlO 6 



5 10 

A„, [km/s] 



5 10 
A rol [km/s] 



Notes. The meaning of the symbols is the same as in Tab. [4] 
" LCl=(fi - R)/(B + R + V) + 0.47M V , is the linear combination of the 
pair of listed parameters displaying the strongest correlation with A rol . 
b hC2=M v - 1.73[Fe/H], displaying a correlation with A rot nearly as 
strong as LC 1 . 



Fig. 11. A ro , vs. the four parameters displaying the strongest cor- 
relation with it (Tab. B). In all cases, removing NGC6441 from 
the sample has a negligible effect on the final value of P,. 

Table 4. Spearman rank correlation coefficients for A ro ,/cro 



Parameter 


Spearman s a 


W 




Pf 


(B-R)/(B+R+V) 


-0.769 


5.526 


23 


2x1 0~ 5 


[Fe/H] 


0.492 


2.652 


24 


0.015 


IQR 


-0.336 


1.595 


22 


0.126 


Z 


0.271 


1.322 


24 


0.200 


Roc 


-0.211 


1.014 


24 


0.321 


logpo 


-0.196 


0.936 


24 


0.359 


7* 

mti.x 


-0.188 


0.835 


21 


0.414 


R *per 

ell 


-0.187 


0.829 


21 


0.417 


0.145 


0.673 


23 


0.508 


My 


-0.076 


0.360 


24 


0.722 


<T 


-0.010 


0.049 


24 


0.961 



Notes. " Spearman rank correlation coefficient of A TOf /cr vs. the con- 
sidered parameter. 

* Absolute value of the associated t variable as defined in Pres s et all 
(1992). 

c Number of clusters involved in the correlation, t is approximately dis- 
tributed as a Student distribution with N c i us —2 degrees of freedom ([Press] 
|etal.|1992| >. 

" Probability to obtain a value of s higher than or equal to observed, in 
absolute value [i.e. P(\s\ > |j' („|)] from two independent variables (i.e. 
in absence of correlation), from a Student test. 

* Peri-galactic distance (R per ) and maximum distance above the Galactic 
plane (Z max ) are taken from |Dinescu et aL| (Tl999 ); Casetti-Dinesc uet al.| 
((2007), as in |Carretta et"aI7|<2010a| l. ~ 
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Fig. A.l. Left panels: probability distribution of observing a ro- 
tation amplitude (A bs) given a true value of the amplitude (A, rue ) 
and a uniform distribution of inclination angles (z). Different 
tones of gray indicate different probability levels, from 8% 
(darkest gray) to 1% (lightest gray), in steps of 1%; white re- 
gions corresponds to probability below 1%. The continuous line 
is the A a bs = A true relation. In the upper panel A rrue values are 
drawn from a uniform distribution, in the lower panel they are 
drawn from a Gaussian distribution. The probability is estimated 
in square bins with sides of 1.0 km s _1 spaced by 0.25 km s _1 in 
both directions of the plane. Right Panel: distribution of sin(z') for 
a uniform distribution of i. 



Appendix A: The projected rotation amplitude as a 
tracer of the true rotation amplitude 

To explore quantitatively the effect of projection on the distri- 
bution of observed rotation amplitudes of a population of GCs 
(or any other system whose inclination w.r.t. the plane of the sky 
z is unknown and unconstrained), we performed the following 
simulation. In a cluster with a true rotation amplitude A true , we 
detect an observed rotation amplitude A„b s = A true sm(i), where 
0° < i < 90° (z = 0° corresponds to a cluster seen pole-on, 
i = 90° corresponds to a cluster seen "edge-on", i.e. having 

Aobs — Af rue ). 

We simulated a population of 10 6 clusters whose inclination 
is drawn from a uniform distribution, since there is no preferred 
value for z. We explored two different distributions of A, rue : a uni- 
form distribution in the range 0.0 km s~'< A true < 10.0 km s , 
and a Gaussian distribution with mean (A, rue ) =4.0 km s'and 
cr =2.0 km s~ 



In the lefthand panels of Fig. A. 1 we plot the 
density distributions obtained in the two cases in the plane A true 
vs. A bs- It is quite clear that the regions close to the A b s = A true 
relation display the highest probability density. Independently 
of the assumed distribution of A true , for a uniform distribution of 
inclinations, observed rotation amplitudes are thus more likely 
close to the true value. 

This is due to the simple fact that the sinus function maps a 
uniform distribution of angles i between 0° and 90° into a dis- 
tribution of sin(z) that is strongly peaked toward sin(z)=1.0, as 
shown in the upper righthand panel of Fig. A.l In particular 



sin(z')> 0.5 in ^66% of the cases, sin(z)> 0.7 in ^50% of the 
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cases, and sin(;)> 0.9 in ^28% of the cases, while sin(;)< 0.2 in 
just ^12% of the cases. We can therefore conclude that A„h s is 
a reasonable proxy for A true in a statistical sense. This supports 
the idea that the observed correlations between A rol or A rot /(To 
and, e.g., My, [Fe/H] or B ^_ v , may trace physical correlations 
between these parameters and the true rotational amplitude. 
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